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Abstract 

A standard high purity germanium (HPGe) detector was used to determine the 
previously unknown neutron activity of a weak americium-beryllium (AmBe) 
neutron source. 7 rays were created through 27 Al(n,n'), 27 Al(n, 7 ) and 1 H(n, 7 ) 
reactions induced by the neutrons on aluminum and acrylic disks, respectively. 
These 7 rays were measured using the HPGe detector. Given the unorthodox 
experimental arrangement, a Monte Carlo simulation was developed to model 
the efficiency of the detector system to determine the neutron activity from 
the measured 7 rays. The activity of our neutron source was determined to be 
307.4 ± 5.0 n/s and is consistent for the different neutron-induced 7 rays. 

Keywords: neutron activation, germanium detector, simulation, spectroscopy, 
activity determination 


1. Introduction 

As neutrons are difficult to detect, determining the absolute activity of a 
neutron source is challenging. This difficulty increases as the activity of the 
source decreases. Sophisticated techniques exist for neutron activity measure¬ 
ments, including the manganese bath technique^, proton recoil techniques^!] 
and the use of 3 He proportional counters[3]; nevertheless, the development of a 
method utilizing commonly available high purity germanium (HPGe) detectors 
would be advantageous. 

HPGe’s are an industry standard for measuring 7 ray energies to high preci¬ 
sion. Neutron Activation Analysis [21 [5] (NAA) is ordinarily used to determine 
the elemental composition of materials, as well as to perform neutron flux and 
activity determinations. When a target material is exposed to a neutron flux, 
it results in the activation of previously stable isotopes in the material. The 
7 rays from the subsequent decay of the unstable activation products may be 
measured and used to determine the neutron flux or activity. 

When determining neutron source activities on the order of 10 2 n/s, as in 
this study, the use of traditional NAA is inadequate^]. The low neutron flux 
produces insufficient quantities of radionuclides, resulting in levels of radioac¬ 
tivity well below the detection limits of standard 7 ray spectrometer systems. 
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Figure 1: Scale representation of the major components of the experimental 
detector system. The inner and outer shielding are shown in a cross sectional 
view. 


This paper presents modifications that enable accurate measurements of the 
activity of such weak neutron sources by placing the neutron source and tar¬ 
gets simultaneously on the HPGe detector combined with use of a Monte Carlo 
simulation to determine the production rate and detection efficiency of 7 rays 
produced from the neutron interactions. 

2. Experimental Setup 

The detector system consists of a GC12023 HPGe detector from Canberra 
Industries. The crystal is 87.1 mm in diameter by 90.1 mm in length and 
is located within a 101.6 mm diameter aluminum end cap. The detector is 
arranged at the centre of a graded rectangular shield. The inner shield consists 
of 25.4 mm thick copper and the outer shield of 250 mm thick lead. These 
components are shown in Figure |T| A Canberra Lynx digital signal analyzer 
acquires the spectral data. The response of the detector is well constrained 
through calibrations using a standard suite of 7 ray calibration sources. 

The neutron source, made by Eckert & Ziegler Isotope Products, Inc., is an 
americium-beryllium (AmBe) AM241SNA02 source. It consists of a homoge¬ 
neous mixture of americium oxide and beryllium metal within a double stainless 
steel encapsulation, each layer with a thickness of 0.8 mm. The manufacturer- 
specified activity of the 241 Am in the AmBe source is 5.5 MBq; however the 
neutron activity is unknown. Neutrons with energies ranging up to 11 MeV are 
produced via the 9 Be(a,n) reaction with 241 Am as the a source [7j- 

Aluminum and acrylic were used as neutron targets. Aluminum has both fast 
( 2 'Al(n,n') and 27 Al(n,p)) and thermal ( 2 'Al(n, 7 )) neutron reactions. Acrylic 
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Figure 2: Cross sections of the various neutron reactions^. The neutron energy 
spectrum of an AmBe source is included for reference [7], 


serves well as a neutron moderator, with its high hydrogen content, and also 
has a thermal neutron reaction ( 1 H(n, 7 )). The cross sections for these reactions 
are shown in Figure [2] together with the energy spectrum of neutrons produced 
from an AmBe source. 

The targets were a 17.7 cm diameter by 2.7 cm thick aluminum disk (1.8 kg) 
and a set of six 19.0 cm diameter by 1.2 cm thick acrylic disks (0.4 kg each), 
with one disk having a 1.0 cm diameter hole in the centre to accommodate 
the neutron source, when necessary. The AmBe source irradiated these targets 
generating the 7 rays from the reactions mentioned. Because of its low activity 
and the short lifetimes of the reaction products, the neutron source was placed 
along with the targets in the HPGe detector system. 

2.1. Measurements 

Each measurement (with the AmBe source, target(s), and HPGe) results in 
a spectrum of gamma ray energies that contain photopeaks at the energies cor¬ 
responding to specific neutron reactions. Each peak has three components: 1) 
7 rays originating from the AmBe source neutron interactions in the target ma¬ 
terials; 2) background 7 rays (either inherent in the detector system and AmBe 
source, or from AmBe source neutron interactions with non-target materials); 
and 3) background neutrons naturally present in the laboratory interacting with 
the target materials. To determine the activity of the AmBe source alone, the 
contributions from the background 7 rays and laboratory neutrons must be 
removed. This requires 3 measurements per geometric configuration. 

In principle, the neutron activity of the AmBe source may be determined 
from one set of 3 measurements, and by analyzing one specific 7 ray resulting 
from one specific neutron reaction. However, we use a total of 7 geometric 
configurations and 3 different neutron reactions to more thoroughly evaluate 
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Figure 3: Experimental setup of geometric configuration 5. A) shows the ac¬ 
tivation and measurement of the aluminum disk. B) shows the corresponding 
background measurement. The number and type of disks were varied. These 
target materials were always placed between the source and detector. 


the method presented. This allows us to later evaluate the consistency in the 
neutron activities determined by this method and determine a more precise 
neutron activity through a combined fit. 

Table [l] lists the geometric configurations for all measurements performed, 
each for a period of approximately 24 hours. Varying the configurations (the 
geometry) in this way enabled the different thermal and fast neutron reactions 
to be measured independently. Each reaction produced distinct peaks in the 
7 ray energy spectrum. The count rate of any particular peak, A react i on , is 
determined by: 

Areaction — Ameas Abkgd Ai a b- ( 1 ) 

A m eas is the peak count rate in a configuration with the target and source, Abkgd, 
with the source only, and Ai a b with the target only. These correspond to columns 
2, 3 and 4 of Table [l] respectively, and each represents a separate measurement 
taken on the HPGe detector. An example diagram of the experimental setup 
used to measure A meas is shown in Figure [3jA. 

Abkgd consists of the 7 ray backgrounds (dominated by those originating 
from the AmBe source) including those resulting from activation of non-target 
components of the detector system, such as the aluminum end cap of the HPGe. 
These measurements were done in a way such that the only volume missing was 
the target volume. In the case of the aluminum disk geometries (configurations 
1-5 in Table [I]), the disk was removed and replaced with thin spacers to elevate 
the remaining acrylic disks to maintain the neutron source geometry, an example 
of which is shown in Figure [3j3. For the acrylic disk geometries (configurations 
6 and 7 in Table [l]), the disks were removed and the neutron source was placed 
on a spacer only. 

Ai a b is the contribution from background neutrons in the laboratory, such 
as those from cosmogenic sources. These rates were determined by performing 
the measurements without the neutron source present. A small contribution 
was observed from the 1 H(n, 7 ) reaction in geometries 6 and 7 (corresponding 
to 1.2%-1.6% of A meas ). The effect of Ai a b was negligible for all of the other 
reactions across the series of geometries. 

2.2. Damage Considerations 

This prompt gamma NAA (PGNAA) approach, with its simultaneous ac¬ 
tivation and measurement, exposes the HPGe crystal to a small neutron flux 
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Table 1: Geometric configurations of the measurements performed. The disks 
and source were placed on top of the HPGe detector in the order listed. Three 
measurements are made for each geometry ID in order to isolate the contribution 
from reactions on the target alone. 


ID 

Measurement 

Background from AmBe 

Background from laboratory 

1 

A1 disk, AmBe 

Spacer, AmBe 

A1 disk 

2 

A1 disk, 1 acrylic disk, AmBe 

Spacer, 1 acrylic disk, AmBe 

A1 disk, 1 acrylic disk 

3 

A1 disk, 2 acrylic disks, AmBe 

Spacer, 2 acrylic disks, AmBe 

A1 disk, 2 acrylic disks 

4 

A1 disk, 3 acrylic disks, AmBe 

Spacer, 3 acrylic disks, AmBe 

A1 disk, 3 acrylic disks 

5 

A1 disk, 4 acrylic disks, AmBe 

Spacer, 4 acrylic disks, AmBe 

A1 disk, 4 acrylic disks 

6 

5 acrylic disks, AmBe 

Spacer, AmBe 

5 acrylic disks 

7 

6 acrylic disks, AmBe 

Spacer, AmBe 

6 acrylic disks 


which may potentially cause damage and result in resolution degradation. The 
total neutron exposure of the HPGe crystal during this experiment can be es¬ 
timated by multiplying the neutron flux through the solid angle subtended by 
the detector by the exposure time in that configuration. Summing over all of 
the measurements performed, this exposure was estimated to be at or below 
10 6 n/cm 2 . Noticeable damage is seen to occur after total neutron exposures 
on the order of 10 9 -10 10 n/cm 2 [91. No noticeable deterioration of the HPGe 
detector’s performance was observed. 

3. Simulation 

In order to determine the number of neutrons required to produce a detection 
in the peak region of the 7 ray energy spectrum, a Monte Carlo simulation was 
used. In the case of simple geometric setups and 7 ray interactions alone, 
standard calibration sources are used to determine this key factor. However, 
this is not possible in the case of complex geometries and with the inclusion of 
neutron interactions. 

Simulations were made using Geant4[T0]. Geant4 tracks particles travel¬ 
ing through matter, models their interactions and allows the user to extract 
information of interest. To simulate an experiment, the user must define the 
geometry, physics processes and data handling during tracking. Critical com¬ 
ponents of the Geant4 simulation are described below. These simulations were 
used to convert the measured reaction rates to neutron source strengths. 

3.1. Geometry 

Each experimental geometric configuration was modeled in detail, including 
the copper and lead shielding, HPGe, aluminum and acrylic plates, and neu¬ 
tron source. The aluminum target was defined with the elemental composition 
of the aluminum alloy used (6061), which includes small amounts of magnesium, 
silicon, iron, copper, chromium, and trace amounts of other elements. The alu¬ 
minum concentration in the alloy was determined by conventional Instrumental 
Neutron Activation Analysis (INAA) at the University of Alberta SLOWPOKE 
Nuclear Reactor Facility and found to be 97.4 ± 0.4 wt. %. 

Neutrons were generated isotropically from positions distributed uniformly 
within the central powder of the AmBe geometry. The energies of the neutrons 
were drawn from the distribution in Figure [2] 
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3.2. Physics 

Within Geant version 4.10.0, the “Standard Electromagnetic” package was 
used to handle the e~ and 7 ray transport and interactions. The “Radioac¬ 
tive Decay” physics list was used to model the decay of the isotopes generated 
through neutron interactions. This includes the generation of all decay prod¬ 
ucts (e“’s, 7 rays, metastable isotopes, etc.) up until a stable isotope state is 
reached. The “QGSP_BIC_HP” hadronic physics list uses the “High Precision” 
neutron cross sections to handle the neutron transport and reactions in the nom¬ 
inal 0-10 MeV energy range. The inelastic scattering process was overridden to 
produce excited aluminum atoms which were then handled by the radioactive 
decay processes to produce 7 rays. This procedure resulted in 7 ray spectra 
consistent in peak width and location with measured spectra. 


3.3. Data Handling 

An event is defined as the generation of one neutron along with all of the 
secondary particles it produces along its trajectory through the simulated ex¬ 
perimental geometry. Several different values are extracted and recorded during 
each of the events. Most importantly, the total energy deposited in the HPGe 
detector (primarily from 7 rays) for each neutron generated is recorded. 

10 7 neutrons ( N S i mu i ate d ) were simulated for each experimental geometry. A 
7 ray spectrum was created by summing over the total energy deposited in each 
event. The same analysis algorithm was used to determine the integrals under 
the peaks of the simulated spectra as was used for the spectra produced from 
measurements. Integrating under the peaks ( N pea k ) in these simulated spectra, 
the neutron source strength (A n ) was determined from the measured data as: 


A n — A r 


N. 


simulated 


N, 


peak 


( 2 ) 


The ratio N pea k/N S i mu i at ed represents the absolute system efficiency, which 
ranged from 1.57xl0 _ 5 -4.94xl0~ 4 , depending on the 7 ray energy. 


4. Calibration 

How accurately the simulation reproduces measured 7 ray spectra was eval¬ 
uated by comparing the simulated efficiency of detecting 7 rays produced by a 
radioactive source to the known efficiency of the HPGe detector for detecting 7 
rays from a calibration source containing the same radionuclide. The most pre¬ 
cise calibration source used for this comparison was a 60 Co point source with 
an activity known to 1.9 % uncertainty. 60 Co emits 7 rays at 1173 keV and 
1332 keV (which may sum to 2505 keV) which are in the energy region of the 7 
rays produced by the target neutron reactions. 

This 60 Co source was also moved radially from the centre of the end cap to 
investigate the change in efficiency over the crystal surface (and to gauge the 
accuracy of the modelled internal geometry of the HPGe detector). In all cases, 
the efficiency determined from the simulation was within the uncertainty of the 
measured efficiency from the HPGe detector (dominated by the uncertainty of 
the calibration source activity). This efficiency comparison is shown in Figure [4] 
As a result, we conclude that the simulation models the total efficiency of mea¬ 
suring 7 rays with this specific detector to a precision of 1.9 %. This is taken 
into account in later neutron activity determinations. 
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Figure 4: Measured and simulated efficiencies for detecting 60 Co through either 
of the 1173 keV photopeak, 1332 keV photopeak, or 2505 keV sum peak. The 
efficiency is shown for varying distances of the calibration source from the centre 
of the detector’s end cap. 


5. Results and Discussion 

Performing the series of measurements described in Table |TJ peaks corre¬ 
sponding to the 7 rays from the target reactions were observed. These will be 
discussed in detail and presented alongside the associated simulated measure¬ 
ments. 

5.1. Fast Neutron Reactions 

The 2 , Al(n,n') reaction, which dominates in the fast neutron regime, was 
observed when the aluminum disk was present. The 844 keV and 1015 keV 
peaks, corresponding to the 1 st and 2 nd energy level transitions of 27 A 1 , are 
shown in Figure [5] and are compared to the simulation. The peak corresponding 
to the 3 rd energy level transition at 2212 keV was also observed in the data at 
a much lower activity. The negative peaks observed in the measured spectrum 
result from the background subtraction, where the removal of material results 
in less attenuation of background peaks. 

The 27 Al(n,p) reaction contributes to the 844 keV and 1015 keV peaks. 
Simulations showed that this contribution is small (2.0%-3.0% for the different 
geometries). This contribution is taken into account in the subsequent analysis. 

From evaluation of the time between interactions in the simulation, the fast 
neutron reactions occur <10 ns after neutron emission. For AmBe sources, a 
4438 keV 7 ray is sometimes emitted along with the neutron. This occurs for 
(57.5 ± 2.8) % of the neutrons produced[Il ,. This 7 ray may also deposit energy 
in the HPGe detector in coincidence with an aluminum 7 ray. Simulating this 
4438 keV 7 ray provided a coincidence correction to the total efficiency applied 
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Figure 5: Measured and simulated energy spectra for the counting in geometric 
configuration 5 (four acrylic disks on top of the large aluminum disk). The 
peaks seen are the result of the 27 Al(n,n') reaction. 


to a given photopeak. This correction for the 27 Al(n,n / ) was small, ranging from 
3.2-6.1 %. 

5.2. Thermal Neutron Reactions 

Peaks associated with the 27 Al(n, 7 ) and 1 H(n, 7 ) reactions were observed 
for experimental geometries 3-7. Figure [ 6 ] illustrates these results and also 
demonstrates the effectiveness of acrylic as a neutron moderator. 

The 1 H(n, 7 ) 2225 keV peak was seen in all spectra from all geometries with 
acrylic. A peak from geometry 7 is shown in Figure [7J The 4438 keV 7 ray 
correction was negligible for thermal neutron reactions. This is due to the fact 
that the approximate time between the 1 H(n, 7 ) and 2 'Al(n, 7 ) reactions and 
neutron emission was ~10-50 p.s (determined through interaction times in the 
simulation). The 7 rays from these reactions are then detected separately in 
the HPGe detector from the 4438 keV 7 ray emitted from the AmBe source. 

5.3. Activity Determination 

By comparing a measurement to its associated simulation, a neutron source 
activity was calculated for each reaction observed in each geometry. For the 
27 Al(n,n') reaction, the interference-free 1015 keV peak was used for activity 
calculations though multiple peaks were observed. From the full system efficien¬ 
cies determined by the simulations, the neutron source activities were calculated 
from Equations [l] and [5] The results are shown in Figure [SJ A best fit to these 
data yields a neutron activity of 307.4 ± 5.0 n/s. 
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Figure 6: Measured 27 Al(n, 7 ) peak for different counting geometries. An in¬ 
creasing peak activity is seen for increased neutron moderation by the acrylic 
(error bars removed for clarity). 
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Figure 7: Measured and simulated energy spectra for geometric configuration 7 
(six acrylic disks). The peak observed is the result of the 1 H(n, 7 ) reaction. 
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Figure 8: Neutron source strengths extracted from the three reactions. The 
dashed line is fit to the data and corresponds to a neutron activity of 
307.4 ± 5.0 n/s. 

Table 2: Effect on the various systematic evaluated. The range corresponds 
to these effects over the applicable geometric configurations which result in the 
reactions listed. 


Systematic 

Effect on neutron activity 

Source/target positioning 2 'Al(n, 7 ) 

5.3-19.1% 

Source/target positioning 1 H(n, 7 ) 

5.9-6.2% 

Source/target positioning 27 Al(n,n') 

1.8-3.3% 

Neutron spectrum 27 Al(n, 7 ) 

0.8-3.0% 

Neutron spectrum 1 H(n, 7 ) 

0.8-2.9% 

Neutron spectrum 27 Al(n,n') 

0.1-3.2% 


5-4■ Systematics 

A further series of measurements were performed to quantify the level of 
systematic uncertainty present in this experiment. The contributions of the 
different major systematics evaluated are quantified in Table [2j One dominant 
systematic uncertainty resulted from the placement of the different targets and 
neutron source with respect to the HPGe (an uncertainty of 5 mm from nominal 
centre was taken in radial position). 

Further evaluation quantified the uncertainty resulting from our knowledge 
of the neutron source spectrum. Small shifts (up to 10%) were added to different 
parts of the source spectrum with simulations being repeated. The subsequent 
effect was seen to be non-negligible but smaller than that of positioning. The 
total uncertainties attributed to these systematic effects are included along with 
the statistical uncertainties in Figure [Sj 
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5.5. Sensitivity 

For each peak measured in all geometries, the minimum detectable activ¬ 
ity (MDA) was determined from the corresponding background measurements. 
From the aluminum reactions, geometry 1 was observed to have the lowest MDA 
(specifically from the 27 Al(n,n') 1015 keV reaction). This activity corresponds 
to a neutron source activity of 9.5 ± 0.4 n/s. Similarly for acrylic, geometry 7 
had the lowest MDA, which corresponds to a neutron activity of 4.2 ± 0.3 n/s. 
These sensitivities are well below the measured neutron activity of the ArnBe 
source. 


5.6. Activity Determination from the 4438 keV j Ray 

The neutron activity was also determined in a separate analysis by directly 
measuring the 4438 keV 7 ray emitted by the neutron source. To obtain a 
photopeak detection efficiency, the same simulation program was used to model 
the detection of these 7 rays. Figure [9] compares the measured and simulated 
spectra. From the peak activity, the neutron activity of the AmBe source was 


determined to be 305.0 ± 16.1 n/s, consistent with the results of Section 5.3 


This relies on knowing the 7 ray to neutron ratio, which we take again to be 
0.575 ± 0.028 dU. _ 

using the 7 rays from the 27 Al(n,n') reaction 


As mentioned in Section 


5.1 


(1014.6 keV) to determine the neutron activity also included the use of this 
ratio. In determining the neutron activity from these 7 rays, the ratio is used 
as a small coincidence correction on the measurement efficiencies. While this 
makes these two analyses not independent, the effect of the correction based on 
the 7 ray to neutron ratio is small (3.2-6.1 %). On the contrary, determining 
the neutron activity from the 4438 keV 7 rays relies entirely on the accuracy 
of this value. The other reactions used to determine the neutron activity are 
independent of this ratio. 

The large uncertainty in the gamma to neutron ratio leads to a less pre¬ 
cise measurement of the neutron source activity than that from the analysis of 
the neutron reactions on the targets. As a result, we quote the latter as our 
measured value and use the former only as a tool for comparison. 


6. Conclusions 

A method was created to determine the previously unknown strength of a 
weak neutron source using a HPGe detector in conjunction with a Monte Carlo 
simulation. Three different reactions in both the thermal and fast neutron 
regimes were used to determine the neutron activity of an AmBe source. The 
independently-determined activities from a series of geometric configurations 
show consistency and result in a neutron activity of 307.4 ± 5.0 n/s. 
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Figure 9: Measured and simulated energy spectra for the AmBe source when 
placed directly on the detector. The 4438 keV peak is seen along with the 
associated single and double escape peaks. 
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